Preconditioned fully implicit PDE solvers for monument conservation
نویسنده
چکیده
Mathematical models for the description, in a quantitative way, of the damages induced on the monuments by the action of specific pollutants are often systems of nonlinear, possibly degenerate, parabolic equations. Although some the asymptotic properties of the solutions are known, for a short window of time, one needs a numerical approximation scheme in order to have a quantitative forecast at any time of interest. In this paper a fully implicit numerical method is proposed, analyzed and numerically tested for parabolic equations of porous media type and on a systems of two PDEs that models the sulfation of marble in monuments. Due to the nonlinear nature of the underlying mathematical model, the use of a fixed point scheme is required and every step implies the solution of large, locally structured, linear systems. A special effort is devoted to the spectral analysis of the relevant matrices and to the design of appropriate iterative or multi-iterative solvers, with special attention to preconditioned Krylov methods and to multigrid procedures. Numerical experiments for the validation of the analysis complement this contribution.
منابع مشابه
Preconditioned implicit PDE solvers for degenerate parabolic equations with applications to monument conservation
We are interested in the mathematical models for the description, in a quantitative way, of the damages induced on the monuments by the action of specific pollutants. A quantitative knowledge of such a degradation is of great importance in precise scheduling of cleaning or deeper restoration works. The analytical study of the solution of the considered models has been conducted and validated by...
متن کاملPreconditioned Implicit Solvers for Nonlinear PDEs in Monument Conservation
Mathematical models for the description, in a quantitative way, of the damages induced on monuments by the action of specific pollutants are often systems of nonlinear, possibly degenerate, parabolic equations. Although some of the asymptotic properties of the solutions are known, one needs a numerical approximation scheme in order to have a quantitative forecast at any time of interest. In thi...
متن کاملTerascale Optimal PDE Simulations (TOPS), An Enabling Technology Center Scientific Discovery Through Advanced Computing: Integrated Software Infrastructure Centers
iii 1 Background and Significance 1 2 Preliminary Studies 3 2.1 PDE Time Integrators . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 2.2 PDE Nonlinear Solvers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 2.3 PDE-constrained Optimization . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 2.4 Linear Solvers . . . ...
متن کاملPreconditioned Conjugate-Residual Solvers for Helmholtz Equations in Nonhydrostatic Models
Numerical integration of the compressible nonhydrostatic equations using semi-implicit techniques is complicated by the need to solve a Helmholtz equation at each time step. The authors present an accurate and efficient technique for solving the Helmholtz equation using a conjugate-residual (CR) algorithm that is accelerated by ADI preconditioners. These preconditioned CR solvers possess four d...
متن کاملFourth-order numerical solution of a fractional PDE with the nonlinear source term in the electroanalytical chemistry
The aim of this paper is to study the high order difference scheme for the solution of a fractional partial differential equation (PDE) in the electroanalytical chemistry. The space fractional derivative is described in the Riemann-Liouville sense. In the proposed scheme we discretize the space derivative with a fourth-order compact scheme and use the Grunwald- Letnikov discretization of the Ri...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2013